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Abstract. The analysis of x-ray reflectivity data from artificial heterostructures 
usually relies on the homogeneity of optical properties of the constituent materials. 
However, when the x-ray energy is tuned to an absorption edge, this homogeneity 
no longer exists. Within the same material, spatial regions containing elements at 
resonance will have optical properties very different from regions without resonating 
sites. In this situation, models assuming homogeneous optical properties throughout 
the material can fail to describe the reflectivity adequately. As we show here, resonant 
soft x-ray reflectivity is sensitive to these variations, even though the wavelength is 
typically large as compared to the atomic distances over which the optical properties 
vary. We have therefore developed a scheme for analyzing resonant soft x-ray 
reflectivity data, which takes the atomic structure of a material into account by 
"slicing" it into atomic planes with characteristic optical properties. Using LaSrMn04 
as an example, we discuss both the theoretical and experimental implications of this 
approach. Our analysis not only allows to determine important structural information 
such as interface terminations and stacking of atomic layers, but also enables to extract 
depth-resolved spectroscopic information with atomic resolution, thus enhancing the 
capability of the technique to study emergent phenomena at surfaces and interfaces. 


1. Introduction 

Specular x-ray reflectivity is one of the work horses for characterizing thin films and 
multilayers. In simple words, the reflectivity is given by interference of x-rays that 
are reflected at the different interfaces realized in such an artificial heterostructure. 
Referring to the reflection of optical light, an interface can be defined as a region in 
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space where there is a change of the refractive index n. Similarly, also in the x-ray 
range even a small change in n will introduce an interface, thus a traveling x-ray 
wave will be reflected. This high interface sensitivity is what allows to accurately 
determine structural properties of heterostructures such as layer thicknesses and 
interface roughnesses by means of x-ray reflectivity. 

Recently, with the advent of synchrotron radiation, the availability of photon 
sources with very high brilliance and tunable energy has opened the frontiers for x-ray 
reflectivity techniques to study additional properties apart from structure. Electronic 
properties, for instance, can be studied by tuning the x-ray photon energies to an 
absorption edge. At these so-called resonant energies, the refractive index depends 
very strongly on the valence shell properties of the resonant scattering centers and 
hence, the sensitivity to spatial variations of the electronic properties is dramatically 
enhanced. This renders resonant x-ray reflectivity (RXR) an ideal tool to study 
electronic properties and phenomena at surfaces and buried interfaces in an element 
specific and non-destructive way. 

The development of RXR was in particular triggered by the recent progress made 
in the atomic scale synthesis of transition metal oxide (TMO) heterostructures. TMOs 
provide perhaps one of the richest and fruitful fields in condensed matter research in 
terms of electronic properties and emerging novel physics [1-4]. Examples of these 
exotic phenomena are, among others, the formation of a two-dimensional electron gas 
at the polar/non-polar interfaces of LaAlOa/SrTiOs [5] or the proximity effects and 
orbital reconstruction in superconductor/ferromagnet (SC/FM) interfaces [ 6 - 8 ]. All 
these properties are closely related to the transition metal (TM) 3 d and oxygen 2 p 
electrons and their interaction with the crystal lattice. 

RXR experiments have therefore in particular been performed at the TM Z/ 2,3 
edges, where the 3 d electrons of the TM are directly probed. In this way, important 
information has been obtained e.g. about the spatial electron density redistribution of 
the Ni 3 d electrons in LaNiOa/LaAlOa multilayers [9] or the Co valence reconstruction 
at a LaCoOs polar him surface [10]. Moreover, employing the x-ray magnetic circular 
dichroism (XMCD) effect, the magnetization prohle of SC/FM interfaces [11], exchanged 
bias systems [12,13] and other multilayers [14—1T| have been studied. 

The analysis and interpretation of reflectivity is commonly done using the Parratt’s 
[18] or the matrix formalisms [19], assuming homogeneous optical properties throughout 
the constituent materials of a heterostructure. Although, this “slab” approach has shown 
to be very successful in describing off-resonant reflectivities, it is not clear if it still holds 
under resonance conditions. This is particularly critical in single crystalline, epitaxial 
TMO-hlms and heterostructures, whose atomic structures typically realize well defined 
lattice planes containing the resonant scatterer. At resonance, these atomic planes will 
interact very differently with the photon beam than the non-resonant regions of the 
material, which immediately raises the question in how far this situation can still be 
described using a single n, i.e., by assuming an optically homogeneous material. These 
effects are particularly important when studying any sort of electronic reconstruction at 
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surfaces and interfaces with RXR, since they are, in fact, expected to occur on atomic 
length scales as well. 

In this report, we investigate in detail the effects in RXR, which are caused by 
the rapid variations of the x-ray optical properties, which are caused by a periodic 
arrangement of the resonating lattice planes along the growth axis. Using a single 
thin film of LaSrMnCR as a practical example, we derive analytical expressions for 
the reflectivity based on the Parratt’s formalism in which the structure of the film 
is considered as layered, i.e., each plane of atoms is considered as a layer. We find 
that variations of n at interatomic distances can have significant effects on the RXR, 
even in the soft x-ray range, where the wavelength of the photons is usually considered 
large as compared to interatomic distances. Indeed, the sensitivity of RXR to the 
atomic structure of a material enables to extract information about a heterostructure 
like interface terminations and stacking sequences, which significantly extends the 
capabilities of RXR. We show that from this approach, spectral information about buried 
interfaces can be extracted and attributed to one or several specific atomic layers. 

2. Slab versus atomic slices: theory 

When calculating the reflectivity, the crystal structure of the film and how it is simplified 
has important consequences for the calculated intensities, especially at resonant 
conditions. Before we start discussing these effects in detail, we first demonstrate how 
significant these effects can be using the model calculations presented in Fig. 1. Here 
we show the theoretical reflectivities for a 8 unit cells (u.c.) thick LaSrMnCU (LSMO) 
film grown on a NdGaCU (NGO) substrate using three different assumptions for the film 
structure. 

In the first model, which will be referred to as “slab” from now on, the conventional 
approach to reflectivity is applied, i.e., the LSMO film and the NGO substrate are 
described as slabs with homogeneous optical/electronic properties given by its refractive 
index n{u)j (cf. Fig.la). In the second type of model, called “atomic slices” in the 
following, additional information from the crystal structure is included. As shown in 
Fig. lb, the Mn02 (MO) and LaSr02 (LSO) atomic layers of LSMO are represented 
by considering them as thin slices with corresponding refractive indices n(co)MO and 

lso and thicknesses of 1.8617 A and 4.5868 A, respectively. The thicknesses of 
these slices were determined using the fractional atom positions in the unit cell [20] and 
the value of the experimental lattice parameter c = 12.897 A of the film obtained from 
x-ray diffraction [21], As we will describe below, the reflectivity is much more sensitive 
to different stacking sequences than to the absolute value of the chosen slice thickness. 
At energies close to the Mn L edges, the refractive index n(u>)j for all the layers in the 
slab as well as in the atomic slices approach were determined using experimental Mn 
scattering factors f and f" determined from x-ray absorption spectra as described in 
the methods section (cf. Sec. 3). Also, the average LSMO density and the total film 
thickness (8 u.c.) is the same for each model calculation. The NGO is again described 
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as homogeneous slab. For the atomic slices description, we further consider two different 
LSO and MO stacking sequences (cf. Fig. 1 c,d) and compare them with the slab model 
in Fig.la. For the sake of simplicity and in order to focus on the differences between 
the models discussed here, roughness will not be considered in the following description. 
Instead we focus on idealized systems with perfectly sharp interfaces. The realistic case 
with finite interface roughness will be discussed in the analysis of our measurement data 
in section 4.2. 

The calculated reflectivities using these three models are shown in Fig. le. As 
one can see, at photon energies away from the Mn L-edge, (cf. curve at 600 eV in 
Fig. le), the reflectivities are similar to one another especially at small q 2 values. In 
strong contrast to the off-resonant region, the calculated reflectivities at resonance differ 
a lot depending on the model used. Furthermore, the atomic slices calculations for the 
different stackings are distinctly different as well (green and blue curves), showing that 
RXR is able to discriminate between different atomic stacking sequences. All reflectivity 
differences between the models become more pronounced at larger q 2 values. 

At first sight the strong sensitivity of RXR in the soft x-ray range to the atomic 
structure of the him is surprising, because the wavelength of soft x-rays is typically 
considered to be large compared to the atomic structure. For this reason the atomic 
layers in perovskite films are usually not taken into account when analyzing soft RXR 
data [9,11,22], In order to better understand the results presented in Fig. le and to 
discuss the discrepancies observed for the above mentioned models, a closer look at the 
atomic slice model is required. In this model we incorporate the atomic structure of 
the him into our analysis using thin homogeneous slices, as illustrated in Fig. 2. In this 
approach, the lattice planes with an area density t] of discrete sites with form factor / 
are approximated by thin homogeneous slices with thickness d ~ 2 A, density p = rj/d 
and refractive index n. 

The refractive index of the latter is given by the optical theorem as n — 1 + 
'^2 i 2iTpirofi/k 2 [23], where the sum is taken over all distinct atoms in the lattice plane. 
Here r 0 and k are the classical electron radius and the x-ray vacuum wave vector, re¬ 
spectively. 

For (j-polarized light with a momentum transfer vector q = (0, q y: q z ) the amplitude 
rehectivity r from a thin slice is in kinematical approximation given by 

r = Qz Qz,s (1 - exp (iq z>a d)) (1) 

q z + Qz, s 

Qz,s = \Jn 2 q 2 - q 2 « q z - cf/q,- (1 ~ n) - q 2 y q?/q 3 z ■ (1 - nf + ... (2) 

Here q z . s is the z component of the momentum transfer vector in the thin slice. At this 
point, we will take only linear orders of (1 — n) into account for the calculation of the 
reflectivity. This is justified because (1 — n) C 1 so that higher order terms only become 
relevant in g 2 S if q z <C |q|. Since we focus here on RXR at large q z , i.e., geometries far 
away from grazing incidence, neglecting higher orders in (1 — n) is justified. For a thin 
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Figure 1. Comparison of calculated reflectivities for an LSMO film using three 
different structural models, (a) homogeneous slab, (b) the LSMO crystal structure 
showing the different LSO and MO atomic slices, (c) and (d) LSMO film modeled as a 
bilayer structure for the LSO/MO and MO/LSO stacking sequences, respectively. N 
refers to the number of bilayers, in this case we have N = 16 = 8 u.c. The average 
density and the total thickness is the same in each case. Note that an inversion of the 
layer structure, i.e, LSO/MO or MO/LSO, has a dramatic impact on the reflectivity. 
The refractive indices of the layers used for the calculation were determined from the 
experimental TEY spectra of a single layer LSMO (cf. section 3). 


slice, we will therefore write down r as 

2 k 2 

r = —(1 - n)(exp(-iq z d/2) - exp(iq z d/2)). (3) 

9z 

The intensity reflectivity is given by |r| 2 . The two phase factors in the above expression 
correspond to the interference of rays reflected from the top and bottom interface at 
±d/2 of the thin slice. This interference term has been introduced artificially by the 
present approximation and does not exist for a single lattice plane. It is therefore 
important to show that this term can be neglected, i.e., to show that the dependency 
of r on d can be neglected. To this end, we express r, using the above expression for n, 

r = -^ r0/ { i -Mi +O | M)1 }, (4 ) 

which shows that the interference effects caused by the two interfaces at ±d/2 do not 
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lattice plane of discrete scatterers: 

positions (xj, y i? z Q ), fomi factor f 



atomic slice: 



centered at z 0 - thickness d. index of refraction n 


Figure 2. Atomic slice model, (a) description of atomic planes of a crystal as thin 
homogeneous slices with finite thickness d. (b) material whose crystal structure consists 
of atomic planes of different composition A and B in the atomic slice representation. 


enter, as long as (q z d ) 2 /24 <C 1, which in turn holds as long as d A/3. This is the 
case for most of the soft RXR measurements where q z < 0.5 A -1 , if d < 5A. A very 
similar, more general result has been obtained in Ref [24]. 

Corresponding to our example LSMO, we now consider a material with two different 
lattice planes separated by d, which are described by two different thin slices A and B, 
respectively (cf. Fig. 2(b)). The total amplitude reflectivity r tot of the whole system 
is the sum of the scattering from all the slices with the corresponding relative phases. 
Using the leading order term of Eq. 4, one obtains for a film with N unit cells 


T’tot 


47rr 0 

Qz 


N-l 


Y.{dAf A + VBf B e^ d } 

u=0 


^2iq,dv 


4vrro f yf 5(yf) \ _ 

q z \ 1 - e^ d 1 + e^ d / v 


(5) 


with A = 2 Nd the thickness of the film, yf = (?7 a/ A + ysf B )/2 describing the average 
scattering strength of the film material and d(yf) = ( QAf A — ?7 b/ B )/2 representing the 
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difference of the scattering strengths of A and B. For q z d <C 1 this can be approximated 
by 


not - 


4Rrr 0 i8(pf) 

Qz l 2 



( 6 ) 


where the term proportional to pf corresponds exactly to what is obtained by describing 
the him as a single homogeneous slab without internal structure. 

But from Eq. 6 it is also clear that the latter description starts to fail as soon as 
S(pf) is not small compared to pf and q z d ~ 1. Indeed, in soft RXR usually one lattice 
plane is at resonance, while the others are not, which means that 5(pf) will be rather 
large. In addition to this, q z max — 0.5A -1 at the transition metal L-edges and d in tran¬ 
sition metal oxides is typically of the order of 2 A, so that at large momentum transfers 
q z d ~ 1. The important result of our analysis is that in the soft x-ray region, the atomic 
structure can affect the RXR signal. In other words, the corresponding phase differences 
can matter and, hence, the wavelength A can not always be considered infinitely large 
as compared to the lattice spacings. Note also that at a certain resonance 5(pf) can be 
large, even for elements with similar atomic number. At resonance conditions it is in 
fact possible to have a finite S(pf) for the same element with different electronic configu¬ 
rations (e.g., valence or spin state) . According to the above discussion the contributions 
to r tot , which originate from the internal atomic structure of the sample, become signif¬ 
icant at resonance and at large momentum transfer. These two regions are exactly the 
most important ones when soft RXR is used for characterizing electronic reconstruction 
phenomena at interfaces. 


Vacuum 
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NdGa0 3 
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Figure 3. Electronic reconstruction of a 7 u.c. LaSrMnCR film, (a) Structural model 
used for the simulation of the reconstruction. For the first unit cell at the film/substrate 
interface Mn (shown in red) has an electronic reconstruction where 50% of Mn is in 
the 3+ state and the rest is reconstructed, (b) Mn scattering factors, red and blue 
show the real and imaginary part of Mn 3+ and of the reconstructed Mn, which was 
constructed to have a Gaussian shape. The real part is obtained from a Kramer-Kronig 
relation. 


To asses the accuracy of the analysis of RXR data in terms of a standard slab 




















model, we generated data sets consisting of reflectivities at energies close to the Mn £ 2,3 
edges for two layer stackings with the same total layer thicknesses of 7 u.c. (cf. figure 
5 a). Additionally, the first LSMO unit cell on top of the NGO substrate was assumed 
to be reconstructed. Figure 3 a shows the atomic slice model with the reconstruction 
for a him with NGO/./Vx(LSO-MO) stacking. For this reconstructed layer, the Mn was 
set to have 50% nominal Mn 3+ scattering factors (red and blue lines in figure 3 b), as 
obtained from experimental XAS (cf. section 3). The other 50% Mn was assumed to 
be reconstructed. The scattering factors for the reconstructed Mn are shown in figure 
3 b, where f"(u ;) was adopted to have a Gaussian line shape with an energy shift with 
respect to the unreconstructed case. /'(a;) is then obtained from a Kramers-Kronig 
relation. As can be observed in Fig. 4 the reconstructed layer has a strong impact on 
the calculated reflectivities, illustrating the high sensitivity of RXR on the atomic scale. 



Energy [eV] Energy [eV] Energy [eV] 


Figure 4. Reflectivity as a function of energy for two selected q 2 vectors showing the 
effect of reconstruction. 

The calculated reflectivities were then fitted with a slab model. The fitting 
parameters were: the overall scaling factor M, the thicknesses of the reconstructed A rec 
and LSMO A lsmo layers and the amount p rec of reconstructed Mn at the interface. The 
results of the fits for the two different stacking orders are shown in figure 5b and c and 
summarized in table 1. As seen in the table, the resulting thicknesses of the different 
layers are close to the starting values within an error of ~ 1-2 A. Also, the information 
regarding the amount of Mn reconstruction is quite close to the original value, however, 
the relative errors for the fitted values of 30-40% are quite significant. 

Based on this analysis, two main conclusions can be drawn. First, fitting RXR 
data using the slab approach can yield a fairly good description of the experiment. 
Notwithstanding, the use of the slab model introduces errors, since it simplifies the real 
material and its reflectivity by overlooking the contribution of the internal structure of 
the him. This can have important consequences: as can be seen in table 1, the fit to 
stacking 1 yields a reasonable y 2 value, indicating a good fit. In contrast to this, the 
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Energy [eV] Energy [eV] Energy [eV] 


Figure 5. Calculated reflectivities from an LSMO film (7 u.c thick) on an NGO 
substrate with an artificial electronic reconstruction at the film/substrate interface. 
For the sake of simplicity, the system is assumed to be perfectly crystalline with sharp 
interfaces. The calculated reflectivities in (a) correspond to NG0/14-[Mn02-LaSr02] 
(stacking 1) and NG0/14-[LaSr02-Mn02] (stacking 2). (b) and (c) show the fitting 
of the calculated reflectivities using a slab approach for stacking 1 and 2, respectively. 
For stacking 1, the slab approach yields a fairly good description of the reflectivity 
at low q z . For stacking 2, the reflectivity lineshape can not be described very well 
with a slab approach. The total fitted dataset consisted of 10 reflectivities at constant 
energy and 10 reflectivities at constant q z , for each polarization. For clarity only few 
reflectivities are shown. 


Parameters 

Start values 

Fit to Stacking 1 

Fit to Stacking 2 

M 

1 

0.81 

1.24 

Alsmo [A] 

77.38 

79.76 

75.18 

A rec [A] 

12.897 

11.53 

13.83 

Prec [-^] 

0.5 

0.45T0.13 

0.54±0.22 

x 2 

- 

0.59 

1.40 


Table 1. Results of the fits of the calculated reflectivities for two different bilayer 
stacking of LSMO using the slab approach. 


resulting y 2 f° r the fit to stacking 2 is about 2.5 times larger than that of stacking 1, 
implying dubious fit results and parameters. The reflectivity curves, calculated in the 
atomic slices approach for stacking 2, therefore cannot be described well by a slab model, 
which neglects the internal structure of the material. Second, by taking into account 
information about the lattice structure and setting up a corresponding atomic slice 
model enables to retrieve important information, such as the stacking sequence, which 
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is lost when the slab approach is implemented. Therefore, a more accurate description 
of the experiment is obtained when utilizing the atomic slices model. 

3. Methods 

A (001) oriented LaSrMnCh (LSMO) him was grown epitaxially on a (110) oriented 
NdGaOs substrate by the pulsed laser deposition technique. Details on the sample 
preparation and characterization can be found elsewhere [ 21 ]. 

X-ray absorption spectroscopy (XAS) and resonant x-ray rehectometry experiments 
have been performed at the 10ID-2 (REIXS) beamline of the Canadian Light 
Source (Saskatoon, Canada) [25], using linearly cr- and 7 r-polarized light. The XAS 
measurements were carried out in the total electron yield (TEY) mode at two different 
scattering geometries and incoming beam polarizations in order to extract absorption 
spectra corresponding to the directions in-plane and out-of-plane, i.e., E || c and E E c 
of the LSMO, respectively. The absorption was measured at the O K, Mn L 2j 3 and 
La M 45 absorption edges. The RXR experiments were carried out at energies around 
the Mn L 2j 3 and La M 45 absorption edges. The reflected intensities were collected in 
the fixed energy (hxE) and the fixed q, (LxQ) modes. The hxE consist of q,-scans 
(6 — 26) carried out at a fixed photon energy, whereas the hxQ refers to energy scans 
at a fixed scattering vector q,. The selected q z vectors correspond to maxima and 
minima of the thickness oscillations taken from the hxE reflectivity curve measured 
at 641 eV. Determination of the imaginary part of the Mn scattering tensor and the 
scalar atomic scattering factors for La, Nd, Ga, Sr and O was done as described in the 
work by Macke et al. [22], This is, the parallel and perpendicular components of the 
absorption spectra obtained from the XAS measurements were scaled to non-resonant 
tabulated values [26]. The X-ray absorption was measured in between the Mn L 2i 3 edge 
and the La M 4j5 edge to avoid any influence of near-edge oscillations on the scaling. The 
real part is then obtained by performing a Kramers-Kronig transformation. Since all 
experiments were performed at room temperature, i.e, well above the Neel temperature 
(in bulk T/v = 127/1 [27]), there is no long-range magnetic ordering of the Mn-moments. 
In this case, the scattering matrix can be assumed to be diagonal [28]. From the atomic 
scattering factors we could calculate the dielectric tensor e and the refractive index 
n(u>) = 1 — <5(a;) + i/3(uj) of the him as shown in figure 6 . The specular reflectivity 
was calculated with the Parratt’s [18] and matrix [19] formalisms using the ReMagX 
suite [29]. 

X-ray photoelectron spectra (XPS) were recorded using PHI VersaProbe 5000 
spectrometer with monochromatic Al Ka (hu = 1486.6 eV) radiation at pass energy of 
23.5 eV. In order to characterize the chemical composition of the him surface, spectra 
were collected as a function of the take oh angle between the surface of the sample and 
the axis of photoelectron detector. Binding energies of the spectra were calibrated with 
an adventitious carbon C Is emission line at 284.8 eV. The XPS spectra were analyzed 
using XPSPEAK 4.1 software after background subtraction by the Tougaard method. 
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Figure 6. The optical constants 5, /3 of LSMO at the Mn 1,2,3 edge. Determination 
was done as described in the work by Macke et al. [22]. 


The shape of the characteristic peaks in all spectra was considered symmetric with a 
combination of 30% Lorentzian-Gaussian profile. 

4. Experimental results and discussion 

In order to determine a realistic model for the studied heterostructure, we characterized 
the film surface experimentally using XPS and determined a realistic parametrization 
for the LSMO/NGO interface. 

4-1. Surface characterization and description of substrate/film interface 

The existence of a SrO layer on the sample surface is revealed by the XPS data shown in 
7a and 7b, where the take off angle dependence for Sr-3d photoemission lines is presented. 
As it can be seen, Sr-3d shows two sets of doublet peaks shifted by ~ 1.1 eV in binding 
energy indicating two different Sr-0 bonds. Peak II shown in blue has been attributed 
to SrO and peak I (yellow) originates from Sr-0 bond in the LSMO structure [30]. Since 
the peak II contribution to the spectrum measured at 20° is larger in comparison with 
the one measured at 45°, and considering that XPS measurements at small angles are 
more surface sensitive, it can be concluded that SrO segregates at the surface. Such 
a SrO segregation layer is commonly found on the surface of manganites [31-35]. We 
therefore included a SrO surface layer into the model for the reflectivity. 

In addition to this segregation layer, a surface adsorption layer on top of the SrO was 
considered based on the fact that the sample has been exposed to air and contaminants 
such as water molecules, carbon etc., can be adsorbed on the sample surface. As a 
simplification, we consider only the scattering from oxygen in the adsorption layer. 

Furthermore, it is known from the growth of Ruddlesden-Popper compounds on 
a substrate that if the c-axis of the thin film is larger than that of the substrate, the 
substrate terraces will lead to antiphase boundaries in the film [36,37]. Such a scenario 
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Figure 7. Angle-resolved XPS spectra of Sr-3d for LaSrMnC>4 thin film deposited on 
NGO (1 1 0) substrate measured at the angle of (a) 20° and (b) 45°. 
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Figure 8. Antiphase boundaries of the LSMO structure when growing on an NGO 
substrate. Once a given thickness is reached, the film grows epitaxially. 


is depicted in figure 8. From the figure it can be seen that it is relatively easy to 
heal antiphase boundaries in LSMO by inserting related compounds like LaMnOa and 
Lai. 5 Sr L 5 Mn 2 07 into the stack. Once the LSMO thickness of the film increases some 
domains become dominating and finally the majority LSMO domain grows epitaxially. 
The presence of a minority phase containing single atomic layers of LaSr0 2 in this 
system was also demonstrated by means of x-ray diffraction. [21] To model such an 
interface, a transition layer was included with a refractive index that was considered 
to be a linear combination of that of NGO and LSMO. For this layer we have defined 
iMnt = (1 — p) ■ IT- ngo + P ' n lsmo, where the factor p was fitted during the analysis of 
the reflectivities. 
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4-2. Fit of the reflectivity data 

The model used for the analysis of the experimental reflectivities consists of an LSMO 
thin film on an NGO substrate with an NGO/LSMO transition layer, a top SrO 
segregation layer, which itself is covered by a contamination layer (cf. section 4.1). 
Figure 10 shows representative experimental reflectivities (red) together with fits (blue) 
corresponding to the slab (Fig. 10a) and atomic slices approach (Fig. 10b and c). The 
fit results are summarized in table 2. During the fitting procedure, all the measured 
reflectivities at both a and 7r polarization were fitted simultaneously in order to get 
a self-consistent result. The fitting parameters were the overall scaling factor M, the 
background intensity as a fraction of the total beam intensity, interface positions zoj 
of all layers, the adsorbate density po at the contamination layer and the “intermixing” 
parameter p of the interface layer (cf. section 4.1). For the slab approach, all the layer 
thicknesses were fit parameters and a Nevot & Croce roughness approximation [38] was 
used at all interfaces. 

In case of the the atomic slices approach, the Nevot & Croce roughness 
approximation can no longer be applied. Instead we use a new approach that allows for 
a smooth transition between different materials, preserves crystallinity and converges to 
Nevot & Croce roughness in the limit of high roughness. This approach is sketched in Fig 
9. Lets assume we have an interface of two layered materials. Material L consists of two 
layers with composition A and B whereas material S consists of layers with composition 
C. The transition from 100% L into 100% S is parametrized using an envelope error 
function f(z ) = erf(z,a) = (<7\/2tt) Jf oo exp(—( 2 /2a 2 )d(. Here a defines the width of 
the transition, i.e., the roughness. We now define the relative abundance of material S 
in a certain position z through the interface as 

S = 1 - erf ((z - zq) d av , a) (7) 

where z is an integer atomic slice number. In our approach the interface location is 
given by zq, where the relative abundance of material L and S is 50%, respectively. d av 
is the average atomic thickness between material L and S. Correspondingly, the relative 
abundance of L is defined as L = 1 — S. Note that this parametrization allows us to 
choose a given film termination at any interface by fixing the zo value. To illustrated the 
interface definition, let us consider the LSMO/SrO interface as an example. The SrO is 
segmented into slices with density pgrO = 5.01 g/crn 3 and thickness dsro = 2.31 A. The 
relative abundance of SrO (S) at a certain position through the interface is calculated 
using Eq. 7. This allows us to determine the abundance of LSMO (L). 

Once L and S are determined one can then calculate new densities and thicknesses, 
i.e., optical constants, for the layers at the interface by summing the contributions of 
the SrO and LSMO layers as follows 

PMn(^) = L- p Mn , Mn0 2 -Layer if z/3 G Z, else />Mn(*) = 0 

d(z) = L■ d M n0 2 + S- dsrO if z/3 G Z 


( 8 ) 

( 9 ) 
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Figure 9. Parameterization of a generic interface according to Eq. 8. Rough interfaces 
of layered compounds can be modeled with the two continuous parameters roughness 
a and termination point zq. 


The definitions for the other atom densities and atom slice thicknesses are set up in the 
same way. During the fitting procedure, the interface positions z 0 j f° r all interfaces j 
can then be fitted thus yielding the interface termination. Correspondingly, the film 
thickness is determined from the relative distance between z 0 j of the top and bottom 
interfaces. 

Like the Nevot & Croce roughness approximation, this method should only be 
applied if the roughness is much smaller than the film thickness, i.e., in the case of sharp 
and well-defined interfaces. We found that this requirement is fulfilled for our system. 
In order to obtain stable fit results, we had to assume the same roughness a for all 
interfaces. Although this obviously needs not to be fulfilled in the real heterostructure, 
the corresponding model fits well to the experimental results (see below), indicating 
that the roughnesses of the different interfaces in the real material are indeed similar. 

Figure 11 shows the resulting y 2 for different terminations ( Zo ) of the LSMO at 
the LSMO/SrO interface. As shown in the figure, the lowest y 2 is obtained when 
the LSMO at the interface with SrO is terminated with an M11O2 atomic layer. In 
comparison, the highest y 2 , i.e., the poorest fit to the data, is obtained when the 
termination is a LaSr02 double layer. The resulting fits for these two extreme cases 
are shown in Fig. 10b and c, respectively and listed in table 2. This analysis allows to 
conclude that the best fit to the experiment is given by the case in which the LSMO 
film is M11O2 terminated. Regarding the interface with the substrate, although the 
NGO/LSMO termination is less sharp there are two main conclusions we get from the 
fit results. First, the thickness of the transition layer is ~9-10 A (cf. table 2), implying 
that the region where antiphase boundaries appear is less than one LSMO unit cell in 
length. This, together with the obtained small roughness value (~2 A), shows that the 
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Figure 10. Selected resonant x-ray reflectivities measured close to the Mn L 2 ,3 edges, 
(a) shows the fit result based on the slab approach, (b) shows a fit result based on 
an atomic slices approach with the lowest y 2 stacking, whereas (c) shows a fit result 
where the LSMO/SrO termination was fixed to be at the LaSr0 2 layer, (d), (e), and (f) 
show the elementary density profile throughout the film thickness for the corresponding 
models. The lines showing the interface between the transition /LSMO and LSMO/SrO 
layers corresponds to the zq value of the error function (Eq. 7) at that interface. 


NGO/LSMO-interface is very sharp . Second, the fits yield that the first layer that 
grows with few antiphase boundaries on the NGO is a single LSO layer. This is better 
seen in the elementary density profiles in Fig. lOe. The lines showing the interface 
between the transition /LSMO and LSMO/SrO layers corresponds to the zo value of the 
error function at that interface. 

Comparing the resulting fits using both slab and atomic slices (Figs. 10 a,d and 
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Parameters 

Slab 

Atomic slices 
Stacking 1 

Atomic slices 
Stacking 2 

M 

2.16 ± 0.05 

2.08 ± 0.05 

1.759 ± 0.032 

Atrans [A] 

8.51 ± 0.13 

9.10 ± 0.25 

9.90 ± 0.21 

Alsmo [A] 

97.92 ± 0.24 

94.64T 0.24 

97.50 ± 0.10 

AsrO [A] 

6.48 ± 0.11 

7.67 ± 0.17 

7.29 ± 0.27 

Ao [A] 

9.46 ± 0.14 

11.09 ± 0.27 

11.83 ± 0.34 

Po [g/cm 3 ] 

1.142 ± 0.025 

1.040 ±0.024 

1.071 ± 0.021 

[A] 

1.47 ± 0.12 

2.11 ± 0.08 

2.36 ± 0.14 

P 

0.427 ± 0.013 

0.548 ± 0.013 

0.359 ± 0.013 

background/1CT 7 

3.12± 0.13 

2.93 ± 0.13 

2.65 ± 0.11 

zq (LSMO/SrO) 

- 

2.81 ± 0.05 

1.7 (hxed) 

x 2 

0.35 

0.33 

0.59 


Table 2. Fit results of experimental reflectivity data with different structural models. 
Note that a model that neglects atomic slices completely is still better than a model 
with wrong stacking. The given errors were taken from the fit covariance matrix. 
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Figure 11. Fit error y 2 for different LSMO terminations zq at the LSMO/SrO 
interface. The lowest y 2 is obtained for a MnCh terminated film. 


b,e, respectively) it is not obvious which yields the best description of the data. From 
a qualitative point of view, both models reproduce most of the features such as the 
thickness oscillations. Also their fit errors y 2 are similar as shown in table 2. Still, 
the fit using atomic layers has important information that is completely lost when the 
slab approach is used. This is the stacking sequence of the LSO and MO bilayers. 
Such information is very import not only for understanding the him growth but also in 
order to explain physical phenomena, which are determined by terminations, such as 
the effects at the LAO/STO interface. This can only be retrieved by describing the him 
in terms atomic layers that form the crystal structure. 
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The power of this approach is better seen by comparing the fits using the atomic 
slices approach with two different stacking sequence of the bilayers, i.e., figures 10b 
and c. Their y 2 values differ considerably, thus showing that only stacking 1 yields a 
proper description of the experiment. So far, our model assumes a single Mn species. 
However, as discussed above, our analysis already shows that there is a transition layer 
between the NGO-substrate and the LSMO-film as well as a SrO-layer at the top of the 
LSMO-film. The Mn-sites in these regions are therefore located in a different chemical 
environment than the ones inside the film . The properties of the Mn valence shell can 
hence be expected to be different in these different regions of the heterostructure. In 
addition to this, also symmetries are broken at interfaces and polar structures like LSMO 
may exhibit a so-called electronic reconstruction. 

To obtain depth-resolved spectroscopic information about the Mn in the various 
layers, we followed a Krarners-Kronig constrained variational fit approach that was 
outlined in great detail by Kuzmenko [39] and subsequently applied to the case of 
resonant X-ray reflectivity [40]. In this approach, local changes of the complex dielectric 
function e are modeled in the imaginary part by a mesh of triangular functions with 
tunable height. The real part of e can then be easily determined by adding the known real 
parts of the fitted triangular functions. Choosing every second data point in energy as a 
mesh, we were able to simultaneously fit the in-plane (e xx = e yy ) and out-of-plane (e zz ) 
component of the dielectric function employing a Gauss-Newton type fit algorithm.The 
computing time needed for the fitting is short and convergence typically occurred within 
a few minutes on a modern desktop computer. 

We searched for changes in the dielectric function at the top and bottom of the 
LSMO layer, as well as at in the transition layer. In each case, the stoichiometric 
parameters which we already introduced earlier and the dielectric function in-plane and 
out-of-plane were fitted simultaneously. We achieved the strongest improvement of the 
fit result when the variational fit approach was applied to the top of the LSMO film and 
when the changes of the dielectric constant were restricted to the topmost M11O2 layer. 
More specifically, the y 2 -value decreased from 0.33 to 0.22, showing the good agreement 
between model and experiment. The result for the Mn L2,3 spectrum of the topmost 
Mn02 layer is shown in Fig. 12. As can be observed in the figure, the lineshape of 
the topmost Mn02-layer is different from the corresponding layer in the film bulk, thus 
implying a change in electronic properties of Mn. This can have multiple reasons. A 
reduction of the Mn valency from Mn 3+ to Mn + ^ 3- ^ for instance may occur as it has 
been reported for similar systems [41]. In addition, changes of chemical environment and 
local breaking of the symmetry are also common factors that induce changes in spectral 
lineshapes. A detailed analysis of the obtained lineshapes requires a detailed theoretical 
analysis, at least in terms of local multiplets models [42], which is beyond the scope 
of the present study and will be addressed in future work. The important result here 
is that our approach indeed yields depth-resolved Mn L2,3 spectra with atomic layer 
resolution. 

Figure 13 shows all the experimental reflectivities for both polarizations, together 
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Energy [eV] Energy [eV] 


Figure 12. The Kramers-Kronig constrained variational fit of the topmost MnC >2 layer 
reveals an electronic reconstruction at this interface. The spectral shape changes are 
caused by a local chemical environment that is different from the film bulk. Comparison 
with [43] indicates that the increased weight at the L 3 edge might be explained by a 
reduction of the Mn valency. 

The reflectivity for this fit is shown in Fig. 13. The variational fit reduced y 2 to 0.22, 
which is an improvement of 32% over the best purely stoichiometric model. 


with the fit obtained with this final model. Although there are some deviations, the 
overall agreement is very good. Quite remarkable is that the polarization dependence 
is nicely captured. 

5. Summary 

We have shown that the atomic structure of a material can influence soft x-ray RXR 
profiles and hence can be very important for the analysis of RXR data. A new ap¬ 
proach for analyzing RXR in the soft x-ray range, which takes into account the atomic 
structure of a material, has been developed and applied to the RXR-analysis of an 
LaSrMn 04 him. The presented modeling in terms of atomic slices not only provides an 
improved description of the experimental data. It also enables to extract important ad¬ 
ditional information like layer termination and stacking sequence of the atomic planes of 
the him. Furthermore, it allows for a reliable extraction of spectral information about 
a specihc layer of atoms inside a thin him. This additional information is essential 
for the understanding of novel electronic phenomena generated at the interfaces and, 
hence, renders RXR an even more powerful experimental tool to investigate artificial 
heterostructures and devices. 
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Figure 13. Final fit result for the measured cr-polarized (top) and 7r-polarized 
(bottom) reflectivities, which includes the variational fit result of the topmost MnC >2 
layer. Energy dependent scans were scaled to allow a good comparison. 
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